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From:  Deborah  A.  Levin 

Subject:  Report  for  AFOSR  Project  Ended  Calender  Year  2002 
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This  memo  is  in  response  to  a  request  from  AFOSR  for  a  final  report  due  at  the  end 

of  2002.  This  material  was  provided  to  the  project  manager  Dr.  Mitat  Birkan  at  his 
contractor’s  meeting  of  August  2002.  It  comprises  the  project  final  report  and  is  presented 
again  here.  Please  forward  it  to  John  McKee  and  any  other  appropriate  people. 

A  main  task  completed  under  this  project  was  to  understand  the  sensitivity  of  chemical 
models  used  in  the  Direct  Simulateion  Monte  Carlo  method  to  the  pair  selection  techniques 
as  well  as  the  use  of  discrete  versus  continuum  energy  levels.  The  results  are  separated  into 
three  groups: 

1.  The  impact  of  the  parameters  in  the  Arrhenius  fit  (A,  B ,  Ea), 

k(T)  =  ATb  exp  -Ea/kT 

on  the  difference  between  the  rate  constant  determined  by  the  original  Arrhenius  fit  and 
that  produced  using  the  TCE  reaction  probabilities  with  the  discrete  internal  energies 
(uncorrected  and  corrected)  was  examined. 

The  impact  was  investigated  for  the  reaction  H20+0— »OH(A)+OH(A),  for  which  A  = 
3.8  x  1CT21  m3/s,  B  =  1.3,  Ea  =  7.69  x  1CT19  J. 

First,  the  parameter  A  does  not  effect  the  difference  between  the  uncorrected  reaction 
rate  and  the  Arrhenius  values. 
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Parameter  B  visibly  impacts  the  difference.  The  more  B  differs  from  0,  the  more  the 
difference.  This  is  shown  in  Fig.  1. 

The  influence  of  the  reaction  threshold,  Ea,  is  even  more  pronounced.  This  is  illustrated 
in  Fig.  2. 

We  can  therefore  say  that  the  uncorrected  rate  leads  to  larger  errors  for  larger  slope  of 
the  reaction  rate  constant  as  function  of  temperature. 

Comparison  of  different  reactants  for  same  reaction  rate  constant  ( A  =  3.8  x  10-21  m3/s, 
B  =  1.3,  Ea  =  7.69  x  10~19  J)  shows  that  for  larger  number  of  internal  degrees  of  freedom 
of  reactants  the  difference  is  generally  somewhat  larger  (even  though  the  impact  is  smaller 
than  for  Ea ).  This  is  shown  in  Fig.  3. 

The  accuracy  of  the  correction  routine  was  checked  for  several  values  of  Ea  and  fixed 
A  and  B  (reaction  H2O+O— >OH(A)+OH(A)).  It  is  seen  that  the  correction  enables  one 
to  quite  accurately  reconstruct  the  reaction  rate  constant  given  by  the  Arrhenius  fit  (see 
Figs.4-6). 

2.  Next  part  is  related  to  the  simulation  of  the  flow  about  a  rocket  forebody  at  50  km  and 
velocity  of  5  km/s.  The  main  goal  here  was  to  see  the  impact  of  the  reaction  rate  correction 
of  the  formation  of  OH.  There  are  two  major  processes  included:  water  dissociation  by  air 
species  and  exchange  reaction  between  water  and  atomic  oxygen. 

General  structure  of  the  flow  is  shown  in  Figs.  7  and  8.  The  important  issue  is  the 
molecular  oxygen  dissociation.  It  is  seen  from  Fig.  9  that  the  molecular  oxygen  dissociates 
strongly  in  the  shock,  and  there  is  much  more  atomic  oxygen  than  O2. 

The  comparison  of  H2O  translational,  vibrational,  and  overall  temperature  profiles  along 
the  stagnation  line  (body  is  at  X  =  0)  is  shown  in  Fig.  10.  The  most  important  parameter 
in  terms  of  reactions  is  the  overall  temperature,  since  the  TCE  model  is  used  here.  The 
overall  temperature  in  the  shock  layer  is  about  5500  K. 

For  this  temperature  regime  (see  also  previous  figure),  the  most  important  reaction  is 
the  exchange  reaction  between  H2O  and  O,  as  it  is  seen  from  the  comparison  of  the  rates  of 
different  OH  formation  mechanisms  shown  in  Fig.  11.  The  exchange  reaction  is  determined 
by  A  =  1.13  x  10f16  m3/s,  B  =  0.  Ea  =  4.818  x  10-2°  J.  For  this  reaction  which  rate  changes 
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very  slowly  with  temperature  the  uncorrected  reaction  rate  is  very  close  to  the 
one.  Therefore,  we  can  expect  results  that  are  very  close  for  the  discrete  con; 
uncorrect.ed  models. 

This  is  exactly  what  was  observed.  The  comparison  of  the  number  Isnsity  | 
H2O  and  OH  along  the  stagnation  line  for  the  discrete  corrected  and  ui: correct: 
and  the  conventional  continuum  model  is  shown  in  Figs  12  and  13. 

3.  The  third  topic  is  modeling  of  two  counter  plumes.  One  consists  of  pure  ) 
the  other  consists  of  pure  O.  The  velocity  of  the  plumes  was  chosen  to  5.5  k  1 
total  relative  speed  is  therefore  11  km/s,  similar  to  the  MirEx  condition::. 

Two  reactions  were  considered, 

H20+0— >0H-f20 

and 

H20+0-^0H(A)+0H 

The  computation  were  performed  separately  for  these  processes  to  exclude  ti 1 
ference. 

Again,  three  models  were  used:  continuum  and  discrete  with  and  without  rea  : 
correction.  The  general  flowfield  structure  is  shown  in  Fig.  14  and  15. 

Comparison  of  the  profiles  of  OH(A’)  (reaction  (1))  and  OH(A)  (reac  ;ion  (2)) 
density  along  the  stagnation  line  shows  that  the  discrete  corrected  model  a  \d  the  c  : 
model  produce  similar  results  (this  is  in  fact  what  was  expected),  wherea  the  un : 
discrete  model  significantly  underpredicts  the  production  rate  of  both  Oh  I  A)  an  I 
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to  original  reaction  rate  constants. 
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Figure  2:  The  impact  of  the  parameter  Ea  in  the  Arrhenius  fit  on  the  ratio  of  the  uncorrected 
to  original  reaction  rate  constants. 
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Figure  3:  The  impact  of  reactants  on  the  ratio  of  the  uncorrected  to  original  reaction  rate 
constants. 
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Figure  4:  Corrected  and  uncorrected  rates  for  Ea  =  1.69  x  10 


Figure  5:  Corrected  and  uncorrected  rates  for  Ea  =  5.69  x  10 
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Figure  6:  Corrected  and  uncorrected  rates  for  Ea  —  9.69  x  10  19  J 
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Figure  7:  Total  number  density  field  about  a  rocket  forebody. 
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Figure  8:  TVanslational  temperature  field  about  a  rocket  forebody. 


Figure  9:  Total  number  density  about  a  rocket  forebody. 
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Figure  10:  Comparison  of  different  H2O  temperatures  along  the  stagnation  line. 


Figure  11:  Reaction  rate  constant  as  a  function  of  temperature  for  different  reactions. 
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Figure  12:  Comparison  of  H20  number  density  profiles  along  the  stagnation  line  for  different 
models. 
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Figure  13:  Comparison  of  OH  number  density  profiles  along  the  stagnation  line  for  different 
models. 
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Figure  15:  Translational  temperature  field  for  the  counter  plume  flow. 
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